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Using a Green's function formulation of the superfiuid current j a , where a momentum q s is 
applied to the Cooper pair, we have calculated j a as a function of q 3 , temperature, and impurity 
scattering for a two-band superconductor. We consider both renormalized BCS and full strong- 
coupling Eliashberg theory. There are two peaks in the current as a function of q s due to the two 
energy scales for the gaps and this can give rise to non-standard behavior for the critical current. 
The critical current j c , which is given as the maximum in j a , can exhibit a kink as a function of 
temperature as the maximum is transferred from one peak to other. Other temperature variations 
are also possible and the universal BCS behavior is violated. The details depend on the material 
parameters of the system, such as the amount of coupling between the bands, the gap anisotropy, 
the Fermi velocities, and the density of states of each band. The Ginzburg-Landau relation between 
jc, the penetration depth Al and thermodynamic critical field H c , is modified. Using Eliashberg 
theory with the electron-phonon spectral densities given from bandstructure calculations, we have 
applied our calculations for j s and j c to the case of MgB2 and find agreement with experiment. 

PACS numbers: 74.20-z,74.70.Ad,74.25.Fy,74.25.Sv 



I. INTRODUCTION 



Two-band superconductivity was first proposed in the 
late 1950's and further studied in the 1960's and 1970's 
as a possible explanation for understanding superconduc- 
tivity in s- and d- band metals. However, it was not un- 
til the discovery of superconductivity at relatively high 
temperature in MgP>2 in 2001, that the most promis- 
ing example of two-band superconductivity was found. 
After nearly four years of effort and the development 
of an extensive literature based on this material, it is 
now evident that this compound is an electron-phonon 
mediated superconductor, but with two distinct bands 
and hence two energy gaps. Indeed, this is firmly estab- 
lished from detailed bandstructure calculations that pro- 
vide the electron-phonon spectral densities which, when 
in turn are used in a two-band Eliashberg formalism, 
give rise to predictions which arc in accord with experi- 
ment. In fact, the picture that emerges is one of excel- 
lent agreement with the data, with the match being as 
good as is found between theory and experiment in con- 
ventional one-band superconductors 0. A detailed re- 
view and summary of the level of agreement found for 
the thermodynamic properties, BCS ratios, penetration 
depth and temperature-dependent energy gaps is found 
m Ref. 0. With such an extensive comparison between 
rigorous Eliashberg theory and experiment completed, it 
is now appropriate to move on to the transport properties 
where less theoretical work has been done. 
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In this paper, we consider the critical current which 
is of fundamental interest in the discussion of transport. 
While most technological applications focus on the bulk 
critical current density, which requires understanding of 
flux lattice pinning, etc., here we are presenting the crit- 
ical current in thin films. This geometry is more relevant 
to understanding issues associated with the fundamental 
superconducting state, rather than focussing on material 
details that lead to vortex pinning, etc., and hence can 
be expected to provide further insight into the nature of 
multiband superconductivity in these materials. We note 
in passing, at this point, that some authors use the term 
depairing current for the critical current although, tech- 
nically, the depairing current is the current at which the 
first Cooper pair is broken and the gap decreases, and the 
critical current is the maximum current possible, which 
is greater than the depairing current. In simple one-band 
s-wave BCS theory the values are almost the same and 
hence there is the tendency to use the terms interchange- 
ably. Here, these two quantities can be quite different 
and so we will use the term critical current to refer to 
the global maximum obtainable in the current density. 

The full temperature-dependence of the critical cur- 
rent j c was first studied in conventional one-band s-wave 
BCS superconductors by Rogers^ and Bardeenpj. Fur- 
ther work following on this was done by Parameter [f| 
for the clean limit, and by Maki0,0 for the dirty limit, 
where the latter author introduced a description of the 
superfiuid current density using the Green's function for- 
malism. A study of the effect of of a current on the 
quasiparticle states density of states of a superconductor 
was done by Fulde0. The most general work was that 
of Kupriyanov and LukichevQ, who calculated the full 
temperature dependence of the critical current in BCS 
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theory for arbitrary impurity scattering, which returned 
the usual Ginzburg-Landau results for the critical cur- 
rent near Tl-|lf)j|. Further considerations led to an ini- 
tial attempt to include strong inelastic scattering[llj fol- 
lowed by a full Eliashberg calculation and study of strong 
coupling effects for conventional electron-phonon super- 
conductors done by ourselves|l2j. Thus, the formalism 
and basic one-band results have been well-established in 
both the weak-coupling BCS limit and strong-coupling 
Eliashberg regime, at all temperatures and for arbitrary 
impurity scattering. 

More recently there has been a renewed theoretical 
effort on the critical current as applied to exotic gaps, 
which are anisotropic, such as d-wave [l3l Il4| and f- 
wavefiH ll!^ . The two-band case is, in some sense, a 
highly anisotropic system and hence the study of the 
critical current in this instance also complements these 
recent works. 

With regard to MgB2, specifically, there has been some 
recent experimental work by Kunchur 16], where the crit- 
ical current has been measured in these systems. On the 
theoretical side, Koshelev and Golubov[l7| have gone be- 
yond Ginzburg-Landau to use the Usadel equations to 
study this system and find interesting effects such as a 
kink in the temperature-dependence of the c-axis criti- 
cal current that reflects the underlying two-band nature. 
This type of feature is observed in other properties such 
as the penetration depth which also shows an inflection 
point in the temperature dependence, however, these ef- 
fects may disappear with increased interband coupling. 
Finally, a more recent work has been done on the topic 
of the non-linear term in the superfluid current in MgB2 
and its effect on the nonlinear microwave response, which 
is relevant for device applications p"s| . 

In our work here, we use thegeneral Green's function 
formulation presented by Maki 7] where the Cooper pairs 
are given a finite center of mass momentum q s . This q s 
leads to a boost in the quasiparticle energy (sometimes 
called a Doppler shift). The superfluid velocity v s = 
q s /m is taken as uniform in the case of discussing currents 
in thin films or wires, unlike the Doppler shift effect in 
the context of superfluid currents in the vortex state, 
where v s varies with spatial position. This latter topic 
has been the subject of a large number of papers in the 
recent literature relating to d-wave superconductors|l9j . 

A very fine physical description of the effect of this 
Doppler shift on the quasiparticle density of states and 
the resulting effect on the superfluid density is presented 
in a paper by Xu et al.|20|. where they consider the non- 
linear Meissner effect in d-wave superconductors from a 
starting point that discusses the basics of the superfluid 
current density. 

Thus, with a body of literature developing on the topic 
of critical currents and superfluid currents, it is of interest 
to elucidate the features in the critical current due to 
two-band superconductivity, where unusual effects can 
occur due to a transfer from one band to the other of the 
dominant contribution to superfluid current and hence 



the critical current. 

Our paper is structured in the following manner. In 
Section II, we introduce the basic theoretical equations 
which we have evaluated. We then examine, in Sec- 
tion III, results for renormalized BCS theory near T c and 
at T = and apply these results to MgB2 . In Section IV, 
we present our Eliashberg results for finite temperature 
and with impurities, evaluating the critical current for 
MgB2 and other model parameters. We also compare 
our results to the data for MgB2. Finally, we form our 
conclusions in Section V. 



II. THEORY 

While various approaches have been used in the past 
for calculating the superfluid current, we will use the 
Green's function method as presented by MakiQ- Thus 
the superfluid current density j s is given in terms of finite 
temperature Green's functions as: 

where e is the electron charge, m, the electron mass, p = 
mv is the electron momentum, and T is the temperature. 
The Green's function is given in terms of the Matsubara 
representation by 

G^^Ap) = \iu{ri) -v-q a + £p 3 + A(n)p 1 CT 2 ] _1 , (2) 

where £ = p 2 /2m — fj,, /i is the chemical potential, and 
the pi and Ui are the Pauli matrices for the particle-hole 
and spin spaces, respectively. The q s is the applied su- 
perfluid momentum and it results in a shift of the quasi- 
particle energies, given to first order in q s by v • q s . Here, 
A(n) and uj(n) are the imaginary axis renormalized su- 
perconducting order parameter and renormalized Mat- 
subara frequencies, respectively, and are defined below. 

The integral over energy £ for the combined Eqns. 
and can be done and the result, generalized to two 
bands, is given in terms of a sum of two partial currents 
j s i and j S 2, for the first and second band, respectively: 

js = jsl + js2 (3) 

with 

3en,i t^, r 1 i(u>i(n) ~ is l z)z 

3si = 7tT > / dz =, 

mv Fi n =-oo J -i yj{u>i{n) - i Sl z) 2 + A 2 (n) 

(4) 

where s, = VFiq s , TH is the electron density and Vf% is 
the Fermi velocity of the i'th band (i = 1,2). The non- 
linear Eliashberg equations for Aj(n) = Z i (n)Aj(n) and 
u)i{n) = Zi(n)uj n have been generalized to two bands and 
with the inclusion of the effect of the q s , they are given 
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II E2: 



as 

Aj (n) = 7rTy^y^ [Ay (m - n) - /4j (w c )6»(w c - | w,, 

m j 

" 1 dz Aj(to) 



In this expression, Sj 



i /Zj 



Ay — 2 



(Qj(m) — isjz) 2 + A 2 (m) 



and 



1 dz 



A,-(n) 



Zj, — 1 + Ay , 



(9) 



(10) 



-i ^ « / (oij(n) — isjz) 2 + A|(n) with w D taken to represent either the Debye frequency 

or some other characteristic energy scale representing the 



(5) 



and 



u>i (n) = w n + 7rTy^y^ Ay (m - n) 

m j 

" dz L)j(m) ~ isjZ 



i 2 



(cDj(m) — isjz) 2 + A 2 (m) 



1 



i 2 ./(^W-^+AKn) 



where t£ 



l/(27TTy) and ir. = 1/(2^") are the or- 
dinary and paramagnetic impurity scattering rates, re- 
spectively, and the electron-phonon kernels afjF(Q,) with 



phonon energy f2 enter through 



Xij (m — n) = 2 



Qo! 2 Pij(0) 
fi 2 + (w n - w m ) 



(7) 



Here, the n indexes the n'th Matsubara frequency w„, 
with uj n = (2n - l)irT, where n = 0, ±1, ±2, • • • . The 
are Coulomb repulsions, with a high energy cutoff w c 
usually taken to be about six to ten times the maximum 
phonon frequency. As written, Eqns. Q-lO are for three 
dimensions with isotropic Fermi surface. For general 
anisotropic Fermi surface, a Fermi surface integral would 
remain, which would need to be done numerically^^]. 
This goes beyond the scope of this paper, however, such 
details are not expected to change the qualitative features 
of our results but rather lead to quantitative changes 
only. 

In the following sections, we present both Eliashberg 
results and renormalized BCS (RBCS) results, the latter 
of which can be quite successful at capturing the essen- 
tial features of the Eliashberg calculations without the 
full numerical complications of the more sophisticated 
calculations even for two-band models 0- To develop 
the RBCS results we use the two-square-well approxi- 
mation or X 6S model written for two-bands (for details 
see Ref. 2j), where now 



A,(T, - f 



m,\uj m | <ujd 



E 



dz [Xij - //*,•] A j (T) 



J (uj m - isjz) 2 + A 



(8) 



phonons in the system. Likewise, the current from Eq. Q 
can be reduced to 



3s 



3erii 

mVFi 



- T E 



dz ■ 



i{u> n — lSiZ)z 



^{w n -is lZ ) 2 + Aj{T) 

(11) 

Equations (JSJ and (|llfl will be used in the next section 
to develop analytic results for comparison with full nu- 
merical calculation. 



(6) 



III. RENORMALIZED BCS RESULTS 

A. RBCS near T c 

In order to evaluate the critical current near T c , we 
need to know the effect of the current on the gap near 
T c . In this limit, both A and q s will be small. After long 
but straightforward algebra, to first order in (1 — t) where 
t = T/T c is the reduced temperature, the gap equation 
of Eq. JSJ) reduces to the form 



1 = XuFi 



Xl2 X 



21 



1 — A22F2 



■F l F 2 + A\G{T c )H l , (12) 



with a second equation obtained from Eq. (|12f) by switch- 
ing indices 1 <-> 2. Here 



A, 



Xjj - (4j 

Zi 



(13) 



and 



f F{T ^ hi[ h±pA + 2 - G (T)v* F 2 ql (14) 



where 



J Fi 



= v Fi /(l + X u + Xij) and G{T) = 
-7C(3)/8(ttT) 2 . For A X (T) -> as T -> T c and with 
no superfluid momentum (q s = 0), we recover the equa- 
tion for the critical temperature (symmetric in l<-2)2]: 

1 = (An + A 22 )g + (A12A21 - AiiA 22 )5 2 , (15) 

with g = 1ii(1.13cl>£)/T c ). Throughout this paper, we will 
find it instructive to consider several simplifying limit- 
ing cases. The completely decoupled case corresponds to 
A12 = A21 = and the assumption that An > A 22 . In 
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this case, the first channel determines the T c , which is 
given by 1.13u>d exp(— 1/Au). As an aside, we note that 
the small separable a 2 aniso trop y model studied exten- 
sively in the older literature 23} maps on to this model 
if An = A(l + a) 2 /2, A 22 = A(l - a) 2 /2 and Ai 2 A 2 i = 
A(l— a 2 )/2, with the anisotropy a 2 assumed small. In this 
case, Eq. 1)150 simplifies greatly since A12A21 — A11A22 so 
that g = l/(An+A 22 ) = l/[A(l+a 2 )] = (l+A)/[A(l+a 2 )] 
which gives T c = 1.13wcexp(-(l + A)/[A(l + a 2 )]). This 
is the well-known result that anisotropy increases the 
value of the critical temperature over its isotropic value 
(a 2 = 0). 

Returning to Eq. JT^ which gives the gap just below 
T c , we need to specify Hi. It is given (after more lengthy 
algebra) as 



H x = An + 



ffA i2 A 



21 



l-g\ 



22 



n 2 \ 2 
9 A 21 



(1-5A 22 ) 3 



(16) 



and an equivalent expression for H2 is obtained from 
Eci. (|16fl by switching indices 1 <-> 2. Solving for A 2 (t), 
we find 



*2/ s (1 — t) 1 

G(T C ) xi 3 xi 



2 2 1 



with 



1 

Xi 



1 

Hi 



An 



Ai 2 A 



12^21 



l-g\ 



22 



and 
1 



1 

Hi 



ffAi 2 A 2 
1-.9A 



22 



2ff- 



*2 
; F1 



g 2 A 22 
1 - gA 22 



(17) 



(18) 



*2 
' J F2 



l-g\ 



22 



(19) 

Likewise, the gap A 2 (t), is formed from the above three 
equations by exchanging 1 <-> 2. When q s = 0, these 
equations properly reduce to those obtained previously 
by Nicol and CarbotteQ- In the limit of decoupled 
bands, 1/xi = 1 and 1/xi — v *fi- Some care is re- 
quired when treating the equivalent equations for I/X2 
and l/x 2 because the combination (1 — <?A 22 ) is replaced 
by (1 — 5A11), which is zero. Nevertheless, we can show 
I/X2 = Vx 2 = 0' which is expected on physical grounds 
as the second band does not contribute near T Cl where 
A 2 (i) = Ofor t -»■ 1, and A 2 (t) = (l-t)8(nT c ) 2 /[7((3)], a 
well-known result. We can also recover other well-known 
results for the small separable anisotropy model defined 
previously. Note first that when the Fermi velocities are 
equal, 1/xi m Eq. lfTT3|l reduces to w|?/xii w hh 1/xi 
given by Eq. {!%))■ Further, 1/xi = (1 + a) 2 (l ~ 5a 2 ) and 
l/Xa = (1 - a) 2 (l - 5a 2 ) which leads to (Eq. jTT^l 



A 2 (t) = (l + a) 2 (l-5a 2 ) 



x 8(?rT c ) 2 2 2 ,0' 
(1 - f) -7C(3T-3^ 



(20) 



and an identical equation for the second gap with the 
first factor of (1 + a) 2 replaced by (1 — a) 2 . For q s = 0, 



Eq. I|20|) reduces to the known result that the average 
gap 



A 2 (t) = (l-5a 2 )(l-t) 8 -^ 



(21) 



with Ai(i) = AoO)(l + a) and A 2 (t) = A (t)(l-a). For 
the isotropic case (equivalent to one band) a 2 = and 
Eq. reduces to 



A 2/A n 8(7rT c ) 2 2 2 , 2 



(22) 



as is well known. 

The expression for the contribution to the current from 
the i'th band is given in renormalizcd BCS by Eq. 1)11)1 . 
Expanding near T — T c for small q s and A gives 



3si(t) = 



7 C(3) ^ ■_>,,,., 
m l t'F l 4(7rT c ) 2 



rA 2 (t)^(? s 



(23) 



The total current is the sum over both bands and can be 
written as 



A\(t) 



:A 2 (<) 



7C(3) 
4(^T C )^ 



-eg s 



(24) 



where m* = mj(l + A^ + Ay). Substituting into Eq. 1)24)1 
the new expression from Eq. 1)17)) (and equivalent for 
1 <-> 2) for A 2 (t), one obtains two terms. The first pro- 
portional to q s and the second to g|? . This gives the total 
current j s as a function of g s . The critical current j c {t) 
is obtained as the maximum of j s (t,q s ) as a function of 
q s . After finding the extremum value of q s , we obtain for 

Jc(t) 



Jc(t) 



8enT r 



3^703) 



(l-t) 



3/2 





' m 


1 




1 " 


3/2 


ni 


1 


n 2 


1 


X 




















m 2 


X2_ 




ml 


xi 


m 2 


X2 



1-1/2 

(25) 



We first note that the standard Ginzburg-Landau (1 — 
i) 3 / 2 temperature dependence of the critical current re- 
mains unmodified in the two-band model. As a first check 
on our expression in Eq. 1)25)1 , we can recover the known 
one-band result when both bands are assumed to be iden- 
tical by taking Ay = A/2, Vpi = Vp2, an d the value of 
m = ri2 = tit/2. Here nx is the total electron density 
per unit volume equal to the sum nj +n 2 . We obtain the 
standard Ginzburg-Landau result); 



jc(t) 



8eriT ( ttT c 



3mw F V V?C(3) 



(l-t) 



3/2 



(26) 



Note that all effective mass renormalization (1+A) factors 
have cancelled in Eq. ()26)l . Applying the separable gap 
anisotropy model to Eq. 1)25)1 leads to a modification of 
Eq. 1)26)1 by a multiplicative factor of (1 — 4a 2 ). This gap 
anisotropy will decrease j c (t) over its value for a = 0. 
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Decoupled bands give the usual expression of Eq. I|26|) 
since the second band makes no contribution near T = 
T c , but now n, m, and vp are those of the first band 
only. In the general case, we need to return to Eq. i|25|) 
which tells us how material parameters affect the slope 
of j c (t) 2/3 which varies linearly with (1 — t) near T c . If 
we normalize the current, as we will in the next section, 
to a slope of negative one near T c , this will mean that 
the current at T = will vary with material parameters 
in contrast to the one-band BCS curve for which it is 
a universal number for the clean limit case that we are 
discussing here. 

A well known result of Ginzburg-Landau theory is that 
near T c , H c (t)/X L (t) = (3nV6/c)j c {t), where H c {t) is 
the thermodynamic critical field and Al it) is the London 
penetration depth, both for the case of zero current. In 
a previous paper [2j, we have given expressions for H c (t) 
and 1/A|(i) as t — > 1. Using these with Eq. we 
arrive at the expression relating the two-band H c , Al 
and j c : 



H c (t) _ 3ttV6 



Az(*)jc(*) 



J. 



(27) 



a 2 anisotropy correction to H c (t) and 1/Az,(i) are both 
of the form (1 — 2a 2 ) which gives a factor of (1 — 4a 2 ) 
in the ratio H c (t) / 'Az,(i). This agrees perfectly with the 
anisotropy factor we derived earlier for j c (t). Thus, in 
this case, the Ginzburg-Landau relation will hold. 

We wish to comment on the use of the Ginzburg- 
Landau relation H c (t) / Ai(i) = (3ir\/6/c)j c (t) for com- 
parison with data at all temperatures, which is done by 
substituting two-fluid model forms for H c {t) and A_l(£). 
The Ginzburg-Landau result should only hold near T c 
and indeed we find that when we form the ratio of 
H c {t)/AL.{t) using the H c (t) and A^(i) from the full 
Eliashberg procedure, the resulting j c {t) curve rapidly 
deviates from the Eliashberg calculated j c (t)- The re- 
sults of using such a phenomenological procedure can 
give curves which are lower or higher than the true re- 
sult, and will not even reproduce the qualitative feature 
of a kink in the j c {t) that we have found in the two-band 
calculation. We conclude that the Ginzburg-Landau pro- 
cedure should only be used near T c as is appropriate to 
its derivation and no reliable use can be made of it phe- 
nomenologically for lower temperatures. 



where c is the velocity of light and J is a complicated 
material-dependent parameter. It is given by 




N*v 



*2 
Fl 



N*v 



Xi 



X2 

where we have found it convenient to change from elec- 
tron density n,, to density of states at the Fermi surface 
Ni(0) = Ni. This modifies the usual Ginzburg-Landau 
relation by a constant factor. For the decoupled case 
where 1/ X 2 = l/x' 2 = 0, 1/x'i = vf v and = 1, J 

becomes equal to 1 . There is no change in the Ginzburg- 
Landau relation and this makes sense since near t = 1 
we are dealing with the single band (the first one). For 
the isotropic gap case, where Ay = A/2 for any (i,j), 
from Eqs. (jT5|) and ljT§|l we obtain 1/xi = 1/Xa = 1 an d 
1/x'i = I/X2 = { v *fi + u F2)/2j the average of the squares 
of the Fermi velocities in the two bands, which gives 



J 



2 1 

F2) 



2(N lV * F \ 



N 2 v* F %) 



(29) 



For identical bands, this expression gives 1, as it must. 
This also holds if either the two Fermi velocities are the 
same or the two density of states match, but not if both 
parameters for one band differ from those of the other 
band. We see that, in general, J can be either greater or 
less than one and the Ginzburg-Landau relation ceases 
to hold. 

For identical bands but with gap anisotropy in the sep- 
arable small anisotropy a 2 model, l/x'i — v f/Xi an d 
l/Xa = v f/X2, with 1/xi = (1 + a) 2 (l - 5a 2 ) and 
l/x 2 = (1 — a ) 2 (l ~ 5a 2 ) which gives J = 1 to order 
a 2 . This result is expected since, as is well known, the 



B. RBCS for T = 

We now turn to low temperature and examine results 
of RBCS at T = 0. At zero temperature, the renormal- 
ized BCS gap equation JSJ takes on the form 

Ax = A 11 A 1 a(A 1 ,s 1 ) + A 12 A 2 g(A 2 ,s 2 ) (30) 
A 2 = A 21 A 1 e(A 1 ,s 1 ) + A 22 A 2 S(A 2 ,s 2 ) (31) 

where the barred quantities are renormalized by the fac- 
tor of Zi of Eq. {THl and 



G(A,s) = In 




These equations need to be solved numerically. Results 
are shown in the four lefthand panels of Fig. ^ The 
top frame is for An = 1, A22 = 0.5, A21 = A i2 = 0.01 
(nearly decoupled bands) and Vp\ = vp2- Because A12 
is so small, the larger gap Ai exhibits the known be- 
havior as a function of q s v F1 /A w of a one-band s-wave 
superconductor^ with a small modification. In one- 
band, the gap would be unchanged until q s v F1 / Ai$ ~ 1.0 
at which point it would begin a sharp drop toward zero 
which is reached at q s v F1 / A10 ~ 1.35. Here, the same 
behavior occurs but there is a very slight drop from the 
q s = value of A at q s v F2 / 'A 2 o — 1.0 due to the coupling 
to the second band. The lower gap behaves differently. 
It retains its q s = value until q s v F2 /A 2 o ~ 1.0 at which 
point it also begins a rapid drop but retains a small fi- 
nite value until the point when Ai = 0. This small tail 
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FIG. 1: Renormalized BCS calculations for T=0 showing, 
on the left, the upper and lower gaps as a function of the 
supcrfluid momentum q s . The top frame is for An = 1.0, 
A22 = 0.5 and A12 = A21 = 0.01, with vf2/vfi = 1-0. In the 
second, the interband coupling A12 = A21 has been changed 
to 0.1. In the third frame, the parameters are the same as 
for the first frame but now vf2/vfi = 2.0. Finally, for the 
bottom frame, the parameters are the same as for the top 
frame but now A22 = 0.7. On the right is shown the total 
superfluid current density j s for each case (solid line) and 
those components for the first band j 3 i (dashed line) and and 
the second band j 3 2 (dotted line). 



is due to the coupling A12 = A21 which guarantees that 
the second gap is nonzero, although it can be very small, 
as long as the first gap is finite. In the second frame 
A12 = A21 has been increased to 0.1. This results in 
a considerable integration of the two gaps. At q s = 
the lower gap ratio A2/A10 has increased from the 0.4 
value of the top frame to 0.57. More significantly for 
the present paper, as q s increases both gaps retain their 
q s = value until q s Vp r> /A2 — 1.0, which is the point at 
which the second gap starts to decrease, but this drop 
is not as fast as in the top frame because of the large 
interband coupling. At the same point as the second gap 
starts to decrease, the first is also reduced but by much 
less. Beyond q s v* F1 j A10 — 1.0 both gaps drop sharply as 
expected and reach zero together at q s v* F1 / 'Am ~ 1.27, 
a value which is smaller than the one band value of 1.35. 
This shift in zero along the horizontal axis is also accom- 
panied by a shift to the left, relative to the decoupled 
case, of the inflection points just described, but these 
shifts are small. These are signatures of the two-band 
nature of our system. In the third frame of Fig. ^ (left- 



hand side), we have retained A12 = A21 = 0.01 but now 
have increased the ratio vp^l^Fx to 2.0 instead of 1.0. 
This has the effect of reducing the region of q s in which 
the second gap remains constant at its q s — value and 
is due to the increase in Vp 2 - F° r the last frame, we 
have kept the parameters as in the first frame but have 
increased A22 to a value of 0.7. This has the effect of 
increasing the q s = value of the second gap and also 
increasing the range over which it stays constant. Addi- 
tionally, the tail beyond q s Vp 2 /A 2 o — 1.0 is not as small 
as it is in the top frame due to less gap anisotropy. 
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FIG. 2: Upper frame: The two gaps for MgB2 parameters, 
calculated at T=0 for in the renormalized BCS theory, shown 
as a function of the superfluid momentum q a . Lower frame: 
The density of states N{ui)/N(Q) as a function of uo is shown 
for several values of q s v* F \l 'A10 = 0.0001 (solid), 0.1 (dashed), 
0.2 (dot-dashed) and 0.3 (dotted). 

In renormalized BCS, Eq. I)ll|l for the partial current 
can be evaluated analytically at zero temperature to ob- 
tain 



J«(0) 



1 + Am + A, 



1 + \ vl + \ 



for 



1 - 



1 - 



A»o 
A l0 



< 1 (34) 
2-1 3/2-1 



for— - > 1 

A l0 



(35) 



Results for the current are presented in the righthand 
side frames of Fig. ^ for the four sets of parameters con- 
sidered on the lefthand side. In all cases, the dotted lines 
give the results for the partial contribution to the current 
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coming from the second band and the dashed for those 
from the first. The sum of the two partial contributions 
is the solid curve. A first feature to note is that for the 
three frames for which A12 = A21 = 0.01 (nearly decou- 
pled), the dashed curves are essentially identical. The 
straight line segment has slope 0.5 in our units and the 
zero is very nearly at q s v F1 /Aio = 1.35 as in the one- 
band case. Note that on the vertical axis, we have used 
the total electronic density of states Nt = Ni (0) + N 2 (0) 
with A/i(0) = N 2 (0). If we had used only N^O) then 
the slope of the straight line segment would have been 1 
which corresponds to what is expected in the one-band 
case. A second point to note is that increasing A12 = A21 
to 0.1 as in the second frame from the top has not had 
a large effect on the partial current from the first band 
except to reduce somewhat the value of q s at which the 
current goes to zero. This is accompanied by a small 
shift to the left in the position of the peak and a small 
reduction in its height. By contrast, the effect of the pa- 
rameter changes are much more significant on the dotted 
curve (the current in the second band). In the top frame, 
the slope of the straight line portion of the dotted curve is 
a little larger than 0.5 because of the renormalization fac- 
tors on Vp 2 relative to v F1 . More significantly the partial 
current starts to deviate from linearity at q s v F2 / A 2 q — 1 
at which point it drops rapidly towards zero but remains 
finite (although small) until the current in the first band 
has reached zero. This also applies to the third frame 
although in that case the current in the second band has 
been increased in value because of the increase in VF2 by 
a factor of two. In all cases the total current (solid curve) 
exhibits two peaks. The peak arising at lower values of q s 
has a significant contribution from both bands while the 
other peak at higher values of q s is due almost entirely 
to the partial current of the first band when A12 = A21 
is very small. On the other hand, when there is more 
integration of the two systems (second frame from top) 
both bands can make a significant contribution. In the 
top three frames, the first peak is lower than the sec- 
ond while in the fourth frame, the opposite holds. The 
height of the first peak can be made higher by increasing 
the value of vf2 as in the third frame compared with the 
first (or equivalently, increasing N 2 /N\ as the density of 
states enters in combination with Up and so will have a 
similar effect, which is redundant to discuss here) or by 
increasing A22 as in the fourth, again compared with the 
first. It can also be increased by increasing the off diago- 
nal coupling as in the second frame. When this is done, 
however, the contribution of the second band to the total 
current in the second peak is also increased. Finally we 
note that the peak heights in the two middle frames are 
only very slightly different and the lower q s peak could 
be made to be the highest through a small tuning of the 
parameters. As the highest peak determines the value of 
j c , this transfer from one peak to the other by varying 
parameters is very interesting and will produce nonstan- 
dard effects at finite temperature. 



C. RBCS at T=0 for MgB 2 

Before leaving our discussion of RBCS at T — 0, it 
is of interest to consider the specific case of MgB2 now 
well-established to be a two-band electron-phonon su- 
perconductor. Microscopic parameters have been cal- 
culated from extensions of bandstructure calculations 
to obtain the electron-phonon spectral densities. The 
Coulomb repulsion parameters that enter the gap equa- 
tions (JSJ) and © are also known. The parameters are 
\ aa = 1.017, A W7r = 0.448, A CT7r = 0.213, K a = 0.155, 
M* ff = 0.210, = 0.172, fj,* an = 0.095, fj* = 0.069, 
with lu c — 750 meV. The ratio of the two density of 
states N- K (Q)/N er (Q) — 1.37 and of the Fermi velocities 
vf-k/vfo- = 1.2. To switch to our notation, (er, 71") = (1,2). 
We start with RBCS results. In the top frame of Fig. 
we show the dependence of Ai and A 2 normalized to A10 
as a function of q s v F1 /AiQ. On comparison with some of 
the results of Fig. ^ we note that MgB2 shows consid- 
erable integration of the two bands and also both gaps 
become zero at q s v* F1 / A\q — 1-2 which is considerably 
reduced from the one-band BCS value. 

An interesting quantity that could be measured is the 
quasiparticle density of states N(oj)/N(0) in the presence 
of a current, which was first examined for one-band s- 
wave superconductors by Fulde|||. Such measurements 
have recently been reported for Pb in Ref. [24| . For each 
band separately 



Nj{0) 

LU - 



Y(w + 5i) S 



A? 



(w - Si) 2 



(36) 



with N(oj) = Ni(oj) + N2(u>). In the bottom frame of 
Fig. |3 we give results for N(u)/Nt(0) as a function of 
u)/Aio for various values of q s as indicated in the figure. 
The solid curve is essentially the zero current limit and 
shows the expected superposition of two BCS quasipar- 
ticle density of states with square root singularities at 
the two gaps. As the current is increased in the linear 
region only (cross-referencing with the upper frame of 
Fig. which shows j s for MgB2), the singular structure 
is smeared and the gap region fills as does the region 
between the two gaps. Gapless behavior appears at a 
value of q s near the first peak in the total current due 
to the density of states in the small band going gapless. 
This is in contrast to what is found in a one-band s-wave 
superconductor^ where the gapless behavior occurs at 
the critical current. 

The superfluid current j s versus q s for the RBCS cal- 
culation with MgB2 parameters is plotted in the upper 
frame of Fig. [3J The line types are defined in the same 
way as for the j s curves in Fig. We note that this 
graph is very similar to that shown in the second frame 
(righthand side) of Fig. ^ where it has been discussed 
in relation to other cases. As in that graph, the peak in 
the total current seen at low q s is lower than the one at 



8 



-, i i i | i i i i r 
MgB 2 (RBCS) 





0.6 

o 

<T 0.4 
E 
> 

2 

(S 0.2 



0.0 

0.0 0.5 1.0 

q s V Fl/ A 10 




q s v F/ A io 

FIG. 3: The current j a as a function of q s calculated for the 
MgB2 parameters. Shown are the component contributions of 
the currents for the first and second bands (dashed and dot- 
ted curves, respectively) and the total current (solid curve). 
The upper frame shows the calculation for T — RBCS and 
the bottom frame shows the curves calculated in Eliashberg 
theory for T/T c = 0.1. 



higher q s , so that this second peak determines the critical 
current. It is of considerable interest to compare these 
renormalized BCS results with those from full Eliash- 
berg calculations based on the oq^FiuS) spectra given in 
Ref. |2l| . Very similar values of the parameters can also 
be found in Ref. j2^. Such results are presented in the 
lower frame of Fig. [31 There are several differences that 
are worthy of note. The full Eliashberg results are at 
T/T c = 0.1 and the RBCS are for T = 0. This dif- 
ference in temperature is small and is not expected to 
lead to any significant differences. Examination of the 
partial currents shows that the dashed curves are almost 
the same in both plots except that strong coupling ef- 
fects have shifted the point of zero current to a slightly 
smaller value of q s . This is accompanied by a shift to 
the left of the position of the peak, a reduction in its 
height and an increase in the rounding, all small effects. 
These differences are seen to be slightly more prominent 
in the dotted curve. This shows that RBCS can be used 
with confidence as a first approximation. However, from 
this point on, we will consider only Eliashberg results 
as we have available to us the electron-phonon spectral 
functions and we wish to consider the full effects of tem- 
perature and impurities. 
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FIG. 4: Upper frame: The temperature evolution of the j s 
versus q s for MgB2 calculated in Eliashberg theory. From 
top to bottom, the temperatures range from T/T c = 0.1 to 
0.9 in steps of 0.1. Lower frame: The critical current plot- 
ted as j^ 3 versus T/T c for MgB2 with varying Vf2/vf\ ratio 
equal to 1.2 (solid curve - corresponding to the set of curves 
shown in the upper frame), 2 (short dashed), 3 (dot-dashed), 
4 (dotted), and 5 (long-dashed). The solid dots are the data 
for MgB2 taken from Kunchur[lrj|. The data appears to be in 
good agreement with the solid curve which corresponds to the 
standard set of parameters given for MgB2 in the literature. 



IV. ELIASHBERG CALCULATIONS FOR 
FINITE T AND IMPURITIES 

At the end of the last section, we showed a compari- 
son between RBCS and Eliashberg theory (Fig. EJ) - Given 
that the electron-phonon spectral functions for MgB2 are 
available, it is appropriate to work with the full theory for 
j c as a function of temperature and impurity scattering. 
We begin with a consideration of temperature effects. 
In the upper frame of Fig. we show the current j s (T) 
normalized by 2eiVrw_FiAio/3 as a function of q s Vp 1 / Aio 
for nine values of temperature from T/T c = 0.1 to 0.9 in 
steps of 0.1. The parameters used are those appropriate 
to MgB2 with the upper curve reproduced from the bot- 
tom frame of Fig. |3| As the temperature is increased, 
the peak with larger value of q s remains the highest and 
therefore determines the critical current given as the solid 
curve in the lower frame. Here, j c has been normalized 
by jo so that the slope at t = 1 is negative one by ar- 

2/3 

rangement for the typical critical current plot of j c , 
which is designed to bring out the Ginzburg-Landau be- 
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havior near T c . We see a smooth (concave downward) 
increase in j c (t) as t is reduced. The value of the crit- 
ical current at t = is ~ 0.61. This is lower than the 
one-band BCS value of 0.72. This reduction could reflect 
the two-band nature of MgB2 as well as strong coupling 
corrections 0, E3 due to our use of the full Eliashberg 
equations (J5J and © to calculate the current given by 
Eq. However, reference to the lower frame of Fig. [3] 
shows that the second band contributes only about 10% 
to the total critical current at T — and this remains true 
at T c . Consequently, this curve is very nearly the single 
band result and deviates from the classical BCS curve al- 
most entirely because of retardation effects. Also shown 
is the data of Kunchurflfij. plotted as solid dots. There 
is good agreement between the data and the solid curve 
which is quite remarkable as the calculation has no free 
parameters, in principle, as they are all given in the liter- 
ature. Including impurity scattering has an effect on the 
temperature dependence of j c (t), as we will discuss later 
in relation to Fig. \7\ Unfortunately, information on the 
individual intraband scattering rates in these films for ir 
and a bands, separately, is not available to us. Including 
some impurity scattering in the 7r band would futher im- 
prove the agreement with experiment. The other critical 
current curves in Fig.0]are for different ratios of vfz/vfi 
as indicated in the figure caption. In addition to allow- 
ing for some uncertainty in the reported Dp values in the 
literature, this also mimics the possibility of currents as- 
sociated with the c-axis where the cr-band is known to 
have a much smaller vp than the 7r-band. Likewise, we 
have not included the details of the full Fermi surface 
average in the expression for the current, however, it is 
known that the cr-band is quasi-two-dimensional as op- 
posed to the more three-dimensional 7r-band. Damn and 
Scalapino^l find significant geometrical corrections as- 
sociated with such averaging which could also accentuate 
the difference between the two bands. This may also be 
approximately captured by our use of different vf ratios. 
All the curves are normalized to have a slope of nega- 
tive one at t = 1. As vp2 is increased, the normalized 
critical current at zero temperature increases and can be- 
come considerably larger than its one-band BCS value. 
There is also a change in the temperature variation of 
the curves, the dashed curve shows a clear kink around 
t ~ 0.25 while the long-dashed curve shows concave up- 
ward behavior at intermediate temperatures. This can be 
understood better when the pattern of behavior shown 
in Fig. [3] is considered. The top frame gives the temper- 
ature evolution for the case Vp^/vpi = 2 and the bottom 
is for a ratio of 4. The notation is as for the top frame 
of Fig. 0] and for the same nine reduced temperatures. 
In the top frame, it is the peak at lower q s which de- 
termines the critical current at t = while at higher t, 
it is the peak with higher q s . This crossover occurs be- 
tween t = 0.2 and t = 0.3 and leads to the kink seen in 
the dashed curve of Fig. 0] (lower frame). If we look at 
the individual contributions from each band (not shown 
here) to the total current, we find that at T = it is 
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FIG. 5: The temperature evolution of the j s versus q„ curves 
for MgB2 calculated in Eliashberg theory. From top to bot- 
tom, the temperatures range from T/T c — 0.1 to 0.9 in steps 
of 0.1. The upper frame is for a vfi/vfi ratio equal to 2 and 
the lower frame is for a ratio of 4. 



the second band which makes the dominant contribution 
because of its large value of vf but near T c this contri- 
bution has decayed sufficiently that it only makes a 30% 
contribution to the total. 

In the lower frame of Fig. the peak at smaller value 
of q s has greatly increased in magnitude over its value in 
the top frame. It has also moved to a lower value of q s as 
expected from Fig. 1. However, now there is no crossover 
from smaller to larger q s peak and hence no kink. Never- 
theless, the dotted curve is very different from a one-band 
BCS variation. We can understand better how this can 
arise by once more looking at the partial contribution to 
the critical current from each band. At T = it is the 
second band which dominates while at T c it accounts only 
for 2/3 of the total with the remaining 1/3 contribution 
from band one. This difference in admixture of the two 
bands with temperature is seen to be enough to cause 
an upward curvature in the dotted curve at intermediate 
temperatures as well as to increase the value of the nor- 
malized critical current at T = to a value considerably 
larger than the BCS value of 0.72. Increasing vf2 further 
as in the long-dashed curve of the lower frame of Fig. 0] 
increases the T = value even further as the peak at 
lower q s value becomes even more prominent. Details of 
the temperature variation of the critical current clearly 
depend on details of the microscopic parameters of the 
two bands involved, such as values of Ay, Ni,VFi- For 
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example, it is possible to shift the kink in j c (t) to higher 
values of reduced temperature and to make it sharper 
by considering other parameter sets. To illustrate this, 
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FIG. 6: Upper frame: The temperature evolution of curves 
for j s versus q s for a Lorentzian spectrum with vf2/vfi = 2, 
An = 1, A22 = 0.8, and A12 = A21 = 0.01. From top to 
bottom, the temperatures range from T/T c — 0.1 to 0.9 in 
steps of 0.1. Lower frame: The critical current as a function 
of temperature corresponding to the same curves in the upper 
frame (solid curve) and also for the cases of vf2/vfi = 1 
(dotted) and 3 (dot-dashed). 

in Fig. we show results for a Lorentzian afjF(u)) spec- 
trum defined in detail in our previous paper 0- In this 
model, An = 1 and we take A22 = 0.8 to decrease the gap 
anisotropy as compared with MgB2. We have also taken 
A12 = A21 = 0.01 (nearly decoupled bands). The top 
frame gives the current j s (t) for the case of vf 2 /vfi = 2, 
for the nine values of reduced temperature previously 
chosen. With this vp ratio, the peak with lower value 
of q s remains highest at small t and the crossover of the 
maximum to the other peak occurs at a reduced tem- 
perature between 0.6 and 0.7. This manifests itself as a 
sharp kink in the solid curve for j c (bottom frame) at this 
same temperature. Also shown in this same frame are re- 
sults for the critical current when vf 2 /vf\ — 1 (dotted) 
and vfz/vfx — 3 (dot-dashed). For the larger value of 
VF2, the normalized current at T = is increased over 
the Vf2/vfi = 2 case as expected and for vp 2 = Vpi it 
is reduced but remains above BCS. In this latter case, 
the two peaks in the current versus q s are closer together 
than they are in the top frame of Fig. and seem to 
almost merge before the crossover from lower to higher 



q s peak occurs, which means that the kink in the dotted 
curve can hardly be seen. It is clear that a rich pattern 
of behavior exists for the current versus q s at different re- 
duced t in a two-band model. This translates into quite 
distinct temperature dependence for j c (t) as compared 
with the one-band BCS canonical behavior. 

Some addition insight into this complicated behavior 
can be obtained from a consideration of simplified but 
analytic BCS results in the limiting case of two decou- 
pled, well-separated bands, i.e., for small values of the 
gap anisotropy u = A20/A10. Under this assumption, 
1/Xi = 1, 1/X2 = 0, 1/xi = v% and 1/ X ' 2 = so that 
the critical current near T c given by formula i|25[l reduces 
to 



Ut) = 



16 e7rT r 



AWi(l-i) 3 / 2 , 



(37) 



and is determined entirely by band one, as previously 
commented upon. As the temperature is lowered, how- 
ever, two possible circumstances can arise. The critical 
current can remain determined by the same band (num- 
ber one) in which case we can write approximately 



J'c(0) 
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mivpi 



^10- 



(38) 



A second possibility is that a crossover to the peak at 
lower q s value occurs and we have approximately 



ic(0)-(^- + ^ eA 
mi m 2 l v F2 



20 



For the first case, 



MO) 
jo 
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0.72 



(39) 



(40) 



which is the one-band BCS result where we have normal- 
ized the critical current in such a way that jc^ 3 has a 
slope of negative one at t = 1. This conforms with what 
we have done in the lower frames of Figs. 01 and [|J] For 
the second case, by assumption, there is a crossover and 



Jc(0) 
jo 



2/3 



= 0.72 



A10 j N 2 v F2 , vfi 

u< 1 

1.76T C \N\vfi vf2 



2/3 



(41) 



This simplified approximate formula shows clearly that 
the anisotropy u on its own reduces the value of 
(jc(0)/jo) 2/ ' 3 below its universal BCS value of 0.72 while 
a large value of N 2 relative to N\ or vf2 relative to Vfi 
increases it and that this reduced dimensionless quan- 
tity is no longer universal. It can depend on the band 
structure parameters as well as on the gap anisotropy. 

Next we consider the effect of impurities on the current. 
It is instructive to begin by expanding Eq. J2J| for j s (t) 
to lowest order in q s . Doing so leads to the well-known 
and physically expected result that for small superfluid 
velocity 



(42) 
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where n S i is the superfhiid density of the i'th band given 
for a strong coupling superconductor by 



2ttT 



),-,.! 



Af(n) 



S?(n)+A?(n)]3/a 



(43) 



In BCS at T = 0, this becomes n\/{\ + An + A12) and 
"2/(1 + A22 + A21) for band 1 and band 2, respectively. 
The slope of j s (t)/e as a function of v s gives n s i + n S 2 
with Eq. (|43|l valid for any temperature and any impurity 
content. The impurities enter only through the Eliash- 
berg gap equations 10 and I©- As the temperature is 
increased, the superfluid density is reduced and so the 
slope of j s (t) seen in Figs. 0] to is correspondingly re- 
duced. Impurities also reduce the superfluid density. In 
BCS, the known result for each of the two bands sepa- 
rately to lowest order, is 



1 



1 + x u + x i:j 1 + 4i 



(44) 



8/, 



with £oi = VFi/(7rAio) and U — vpiTu, the coherence 
length and impurity mean free path in the «'th band, re- 
spectively. This applies to intraband scattering which re- 
duces the initial slope of the partial current coming from 
each band. Results are presented in Fig. [7] where we also 
consider the interband case 7^ 0. In the top frame, we 
show j s versus q s for the reduced temperature i = 0.1 
in the Lorentzian model used in Fig. The solid curve 
applies to the pure case and is for comparison. We show 
only three cases, one where tJi = 2T c o (dotted curve), 



= Ci = 0.02T c0 , where 



*22 = ^cO (dashed curve), and t\ 2 = t\ x 
T c q is the T c without impurities (in the case of and 
t 22 , T c will be unaffected, but not so for interband im- 
purities). Intraband impurities reduce only the partial 
contribution to the current coming from that band. It 
reduces its value on the vertical axis and extends the x- 
axis to higher values of q s leaving the second band largely 
unaffected. Of course, since both bands contribute sig- 
nificantly to the total current in the region of the peak 
corresponding to the smaller value of q s , its magnitude 
is reduced even in the dotted curve. On the other hand, 
for the dashed curve, there is little change in the region 
of the second peak which is dominated by the first band. 
By contrast interband impurities modify both bands in 
a similar fashion. The corresponding critical current as 
a function of T/T c is shown in the lower frame. The 
solid line is the pure case of Fig. |B| As expected, a fi- 
nite value of has reduced the critical current below 
its value in the pure case at all temperatures, however, 
when normalized to a slope of magnitude one near T c , 
the effect pushes the curve up at low T . A kink remains 
and the crossover temperature is only slightly shifted to- 
ward higher values. For the dashed curve, the significant 
modifications arise only below the kink temperature. By 
contrast the main effect of the interband impurities is to 
wash out the kink structure leaving the rest of the curve 
largely unaffected (of course the T c would show a small 
shift downwards relative to T c q, not shown here). 
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FIG. 7: Upper frame: The superfluid current j s versus q a for 
a Lorentzian spectrum as used in Fig. [(J where the solid curve 
is repeated from Fig.HJfor T/T c = 0.1 and has a ratio of the 
two vf's equal to 2. The dotted curve is where a scattering 
rate of t n = 2T c o is now introduced, the short dashed is 
for = T c o, and the long-dashed is for tf 2 = £21 ~ 0.02T c o. 
Lower frame: The critical current as a function of temperature 
corresponding to the parameters for the curves shown in the 
upper frame. 



V. CONCLUSIONS 

In conventional one-band s-wave superconductors, the 
energy gap A remains unaffected by the flow of a super- 
current up to a critical value of the superfluid velocity 
Vs = Vs/m (ie., v c s = A/tiivf) at which point pairbreak- 
ing becomes possible. In this region, the current j s is 
linear in v s with proportionality constant equal to the 
superfluid density. Beyond this critical value of v s , the 
rapid reduction in the gap in addition to a quasiparti- 
cle backflow current leads to j s = at A = 0. For two 
nearly decoupled bands, we would expect that the total 
current to exhibit a two-peak structure as a function of 
v s . A first one near the critical value of the superfluid 
velocity for band two and the other near the larger crit- 
ical value of v s corresponding to the first band alone. 
The first peak has in general a significant contribution 
from both bands (unless VF2 3* Vfi) while the second 
peak could be due mainly to the first band. The relative 
admixture depends on microscopic parameters such as 
density of states AT,(0) and Fermi velocities Vpi, as well 
as gap anisotropy. When vp2 is made large as compared 
to Vfi, t ne first peak moves to lower values of v s and can 
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also become larger than the second at low temperature. 
As the temperature is increased the superfluid density 
in each band is reduced but that in the second band is 
more strongly affected because of its smaller gap value. 
Consequently, since the critical current is determined by 
the height of the highest peak, one can have a crossover 
from the first to second peak as T is raised from to T c 
and this will reflect itself as a kink in the temperature 
dependence of the critical current. Such a feature is also 
found in the dirty limit calculations presented in Ref. 
for the critical current along the c-axis. When there is a 
significant interband coupling or the two gaps are close in 
value (isotropic case) the two peaks in the total current 
cease to be independent and a rich pattern of behavior 
emerges for the temperature and impurity dependence 
of the two peaks. In particular, as the temperature is 
raised, the first peak can remain the determining peak 
for j c up to T c or it can become the smaller of the two 
at a crossover temperature which can be increased or de- 
creased through judicious choices of parameters. Another 
possibility is that the two peaks merge into one or that 
the second peak is the dominant one at all temperatures. 
Application to the specific case MgB2, which involves no 
adjustable parameters, gives good agreement with the 
in-plane data of Kunchur [lfij. 

2/3 ■ 

If the resulting critical current j c is normalized to 
have slope of negative one as a function of (1 — t) near 
t = 1, the corresponding value of the critical current at 
T = will differ from its classic one-band BCS value of 
0.72. Strong coupling effects resulting from our use of 
the full Eliashberg equations, with appropriate electron- 
phonon spectral densities describing the elect ron-phonon 
interaction, are known [12j to always reduce the normal- 
ized value of j c (0) while two-band effects can reduce or 
increase it depending on the detailed values of the micro- 
scopic parameters involved. In particular, a large value 
of vf2 relative to vfi can lead to values that can be even 



larger than one. This is expected for the case of current 
orientated along the c-axis for the specific parameters 
representative of MgB 2 . In this case, because the a band 
(band one in our notation) is nearly two-dimensional, 
the Fermi velocity of the 7r band in the c-direction will 
be much greater than that for the a band. 

The superfluid density in each subband, which deter- 
mines the initial slope of the linear in v s regime for the 
supercurrent, is reduced with the increase in intraband 
impurity scattering. This fact can be used to manipulate 
the position along the v s axis and size of the two peaks 
in the total current and consequently the temperature 
dependence of the resulting critical current. Introducing 
impurities in the first band will reduce the critical current 
at all temperatures since this band always makes a con- 
tribution to the total current. On the other hand, if they 
are introduced only in the second band, this will leave 
the critical current largely unaltered around T = T c . In- 
terband impurity scattering affects both bands although 
not necessarily in exactly the same way and leads to the 
smoothing out of the kink in jV^{t) versus t when it 
exists. 

In conclusion, we predict a complex pattern of behavior 
for the total current as a function of v s , temperature and 
impurity content which can be used to restrict further 
the values of the microscopic parameters involved in two- 
band superconductors, such as MgB2. 
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